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o3 ■ Abstract 

S" 

We examine quantitatively the temperature dependence of the linewidths and line shifts in elec- 
tron paramagnetic resonance experiments on single crystals of the single-molecule magnets Fes 



and Mni2, at fixed frequency, with an applied magnetic field along the easy axis. We include 
inter-molecular spin-spin interactions (dipolar and exchange) and distributions in both the uniax- 
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ial anisotropy parameter D and the Lande ^-factor. The temperature dependence of the linewidths 
and the line shifts are mainly caused by the spin-spin interactions. For Fes an d Mni2, the tem- 
perature dependence of the calculated line shifts and linewidths agrees well with the trends of the 
experimental data. The linewidths for Fes reveal a stronger temperature dependence than those for 
Mni2, because for Mni2 a much wider distribution in D overshadows the temperature dependence 
of the spin-spin interactions. For Fes, the line-shift analysis suggests two competing interactions: a 
weak ferromagnetic exchange coupling between neighboring molecules and a longer-ranged dipolar 
interaction. This result could have implications for ordering in Fes at l° w temperatures. 

PACS numbers: 75.50.Xx, 76.30.-v, 75.45. +j 
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I. INTRODUCTION 



Single-molecule magnets (SMM's) consist of identical molecules, each of which made up 
of several magnetic ions surrounded by many different species of atoms. A single molecule 
of the SMM's Mn 12 -acetat©i and Feg^ has an effective ground-state spin of S = 10 and a 
strong crystal-field anisotropy. A zero-field energy barrier against magnetization reversal 
is approximately 65 K (30 K) for uniaxial Mn 12 (biaxial Feg) . 3 i 4 i 5 i 6 i 7 Despite their large ef- 
fective spin, these single-molecule magnets have shown quantum coherence^ and quantum 
tunneling between the energy levels of the two potential wells.— i 10 ' 11 Although dipolar inter- 
actions between different molecules are weak, in the low-temperature limit and near zero 
applied field, a dipolar interaction could stimulate the quantum tunneling and thus explain 
the non-exponential magnetization relaxation observed at early times . 1211 ? 11411 ^ 1 ^ 1 ? 

Recently, electron paramagnetic resonance (EPR) experiments on single crystals of Fes 
and Mn 12 have revealed interesting effects in the widths and positions of the EPR peaks 
as functions of energy level, resonance frequency, and temperature when the applied field 
is along the easy axis . 4 i 18 i 19 i 20 i 21 For both Fe§ and Mn 12 , at fixed frequency, the linewidths 
increase with decreasing energy levels (the largest linewidth corresponding to the transition 
between the ground state and the first excited state), and for a particular transition the 
linewidths increase with decreasing frequency. On the other hand, the details of the tem- 
perature dependence of the linewidths for the two materials are quite different. For Mn 12 , 
the linewidths increase smoothly with increasing temperature, showing a rather weak tem- 
perature dependence. For Feg, for transitions at low resonant fields, the linewidths increase 
sharply with temperature at low temperatures, reach a maximum, and then decrease slowly 
with temperature at higher temperatures. The exception is the transition associated with 
the ground state, for which the linewidth decreases with increasing temperature in the whole 
temperature range studied (2—50 K). On the other hand, for the transitions with high reso- 
nant fields, the linewidths increase monotonically with increasing temperature. For Fes, the 
line positions change non-monotonically with energy level and temperature. 

In Refl20l which was our first attempt to understand the energy level and resonance 
frequency dependence of the linewidths, we showed that for both single crystals of Fes and 
Mni 2 the distribution in the uniaxial anisotropy parameter D of the single-spin Hamil- 
tonian, caused by defects in the samples, contributes substantially to the inhomogeneous 
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linewidths at constant temperature. This was also recently supported by terahertz spec- 
troscopy for Mni2.— The microscopic origin of the distribution in D has not yet been fully 
understood . 23 ! 24 The analysis further showed that for Fe§ the dipolar interactions between 
molecules contribute to the linewidths as significantly as the distribution in D, while for 
Mn 12 their contribution is less significant.— In a recent millimeter- wave study on Feg,— the 
dipolar field may have been overestimated because the distribution in D was not included. 

Since the approximations made in Ref I20I reasonably well explained the linewidth behavior 
at fixed temperature, in the present paper we investigate the temperature dependence of the 
linewidths and line positions. To explain this additional feature, we take into account the 
inter-molecular spin-spin interactions (exchange and dipolar), as well as distributions in the 
uniaxial anisotropy parameter D and the Lande g-factor. We find that the distributions in 
D and g do not contribute to shifts in line positions with temperature, but the exchange and 
dipolar interactions produce nonzero local fields which are temperature dependent (at low 
temperatures), so that they can make the line positions change with temperature. Thus, 
to explain the temperature dependence of the measured line shifts, both the exchange and 
dipolar interactions are needed. Without the exchange interaction, we cannot explain the 
observed non-monotonic line shifts with temperature. This was not included in our earlier 
study, Refl20l. Although the exchange interaction was not considered in our earlier study, 
its effect on the linewidths is negligible as long as the magnitude of the exchange constant 
is much smaller than the linewidths.— 

This paper is organized as follows. A brief summary of the experimental procedures 
is presented in Sec. II. The models for Fes and Mni2 are described separately in Sec. III. 
Sec. IV describes our calculated linewidths and line shifts vs temperature, and their discussed 
in comparison with the experimental data. Our conclusions are provided in Sec. V. 



II. SUMMARY OF EXPERIMENTAL PROCEDURES 

All of the EPR experiments were performed on single crystals with the magnetic field 
aligned along or close to the direction of the easy magnetization axis, in a temperature 
range from 2 K to 50 K. For Mni2, this direction coincides with that of the longest di- 
mension of the essentially needle-shaped crystals. For Fe§, the direction was determined 
by measurement of the largest spread of the resonant field, by rotating the crystal around 
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the approximately known orientation of the easy axis. The single crystals were prepared as 
described earlier . 4 ! 5 ' 27 ' 28 ' 29 EPR measurements were made in the 100—190 GHz range with 
a resonant microwave cavity system described by Hill et al.— , which enables observation of 
distortion-free EPR line shapes. The linewidths were determined by computer-fitting of the 
observed experimental spectra to either a Gaussian or a Lorentzian function for the spectra 
obtained at various temperatures. All the spectra were obtained by keeping the frequency 
fixed and sweeping the field to obtain the resonance peaks, as is usual in EPR spectroscopy.— 



III. MODEL 

For the examined single- crystal sample of Fes,— there are two sources of the EPR line 
shifts: temperature dependence of D and the electronic spin-spin interactions (dipolar and 
exchange interactions) between different molecules. There are also two sources of the EPR 
line broadening: the D-strain effect (distribution in D) and the spin-spin interactions? 2 ^ 
For the examined Mni2 sample,— the sources of the line broadening are the D-strain, the 
^-strain,— and the spin-spin interactions.— 

In our model, hyperfine interactions are not considered for the following reasons: (i) The 
observed line shifts and changes in the linewidths were much larger for Fe§ than for Mn 12 . 
(ii) On the other hand, hyperfine fields are orders of magnitude larger for Mn 12 (in which all 
nuclei have spins of / = 5/2) than for Fe§ (98% of Fe nuclei with 1 = 0). (iii) Any residual 
effect of the hyperfine fields was included in background linewidths. 

Possible mis-alignment of the external field and a spread in the in-plane fields are also 
not considered since mis-alignment of the external field cannot provide the temperature 
dependence of the linewidths and shifts. We checked that for Fe§ the effect of a distribution 
in the transverse anisotropy parameter E on the linewidths was negligible. We also found 
that a spread in the in-plane fields (possibly caused by nuclear spins) can give rise to a 
distribution in the easy axis of each molecule, leading to asymmetries in EPR lineshapes .?? 1 ? 4 



The model and technique used in this study are similar to those in Ref|20|, except that 
we here take into account the non-spherical sample shapes and the actual crystal structures. 
Therefore, we here summarize them only briefly, focusing on what causes the temperature 
dependence of the line shifts and linewidths for each source. The temperature dependence 
of D is discussed in Sec. IV A.l. In our convention, the energy level M s = +10 is the ground 
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state when the field is applied along the positive z axis, while in Refs 
is M s = —10. For clarity, we discuss Fes and Mni 2 separately. 
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the ground state 



A. Fe 8 

We consider an effective single-spin Hamiltonian which satisfies approximate D2 symme- 
try, 

Ho = -DS 2 Z - E(S 2 X - S 2 y ) - g^ B H z S z , (1) 

where the uniaxial anisotropy parameter D = 0.288/c#, the transverse anisotropy parameter 
E = 0.043/cs, 19 g is the Lande g-factor which is close to 2, and //# is the Bohr magneton. 
Here S a is the ath component of the spin angular momentum operator, and H z is the 
longitudinal static applied magnetic field. We assume that the longitudinal magnetic field 
is applied along the easy axis (the z axis), and we ignore the small transverse anisotropy 
terms (E terms) in calculating the linewidths. Thus, the energy level M s is a good quantum 
number of the spin operator S z . According to our convention, the ground state is M s = +10. 

Using the density- matrix equation 35 with the Hamiltonian, Eq. (fl|). and an interaction 
between the spin system and an oscillating transverse field, we calculate the power absorp- 
tion between the levels M s and M s — 1 for a fixed value of the uniaxial anisotropy parameter 
D. In the power absorption, the line-shape function includes a natural linewidth, which is 
a function of temperature. Next we calculate the average power absorption with a Gaus- 
sian distribution in D, where D itself is assumed to be temperature independent. As a 
consequence, the line broadening due to the D-strain effect becomes weakly temperature 
dependent because of the temperature dependence of the natural linewidths. To calculate 
the natural linewidths, we use the strength of the coupling between the spin system and a 



surrounding phonon heat bath obtained in Refl36l For example, for M s = +10, the order 
of magnitude of the natural linewidths is several to several tens of gauss at temperatures 
below several tens of kelvin, and the widths increase with decreasing M s . 

The spin-spin interactions are calculated separately with D fixed, and then combined with 
the D-strain effect to obtain the total linewidths. At low temperatures, all energy levels are 
not equally populated, and the populations of the levels change with temperature according 
to the Boltzmann factor. Thus, the local field on a particular molecule caused by surround- 
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ing molecules changes with temperature. Therefore, the spin-spin interactions are mainly 
responsible for the temperature dependence of both the line shifts and the linewidths . 26 ! 37 To 
calculate the line shifts and the line broadening due to the spin-spin interactions, we use a 
multi-spin Hamiltonian which commutes wit h y^ - 5*|, where the sum runs over all molecules. 
Details of the technique can be found in Refl37. 



n tot = YllHoi + vS)] (2) 



n dipole = W 'Ajufi ■ S k - 3SjSi) , A jk = ^#(3^ - 1) , (4) 
z jk v^vr/ zr jk 

jk 

where 7io« is the single-spin Hamiltonian for the ith molecule, the sum Y^i runs over all 
molecules, and Vi(t) is the interaction between the ith molecule and the oscillating transverse 
field. Here '}i dl P ole is the dipolar interaction between the molecules, and Q k are the direction 
cosines of the vector between molecules j and k (rjk) relative to the easy axis (z axis). 
The sum J2'jk runs over all molecules, so that any two indices are not the same. ~}{ exch is 
the isotropic exchange interaction between the spins of nearest-neighbor molecules where 
the exchange coupling constant is J if the ith and jth spins are nearest neighbors and 
zero otherwise. It is reasonable to assume that FJj V^(i) is much smaller than the dipolar 
and exchange interactions, which again are much smaller than the sum of the single-spin 
Hamiltonians. Since the field is swept at constant frequency, the energy levels change with 
the sweeping field. We neglect slight changes of the energy levels during resonances and use 
H res = (hv — D(2M S — l))/g/iB, where v is an EPR frequency, as the field H in the spin 
Hamiltonian to calculate the energy levels for the particular resonance. 

To calculate the £th moment of the resonant field deviation, we formulate the £th moment 
for a frequency sweep and then convert it to a field sweep. This is justifiable because we 
neglect the slight energy change caused by the change of the field during a resonance. The 
probability density function of the EPR frequency v is given by 



EnE^{exp(-£ n /fc B T) - exp(-£ n ///c B T)} \(n 


J2j Sj 


n')| 2 


J2nJ2n'{ ex P(-^n/k B T) - exp(-£ n , / k B T) } \(n 




n')| 2 



where £ n is the energy eigenvalue of J2j T~(-oj + Ti.^\ \n) is the corresponding eigenvector, yj+, 
denotes the sum over all states \n') such that £ n i > £ n , and J2n J2n' denotes the sum over 



all states for which hu < S n > — £ n < h[y + dv). Using Eq. (jUJ), we calculate the ith. moment, 
, a EnEt>(£n> - £nY{eM-£n/k B T) - eM-£n>/k B T)} \(n\Ej Sf\n')\ 2 



where the temperature dependence of the linewidths is included through the Boltzmann 
factors, and the eigenvalues of TL^ can contribute to (£ n i — En) 1 and/or the Boltzmann 
factors. Assuming that the dipolar and exchange interactions, liP~\ are much smaller than 
the thermal energy, we expand the Boltzmann factor into 



exp 



k B T 



cxp 



E? % 



Oj 



k B T 



k B T 



(8) 



and consider only the first term on the right-hand side. In Ref.38, the Boltzmann factor 
was not included since the temperature of interest was quite high, so that all energy levels 
were equally populated. In our calculations, we use a mean-field approximation, so that the 
sums of Aj k and Jj k [Eqs. (jH) and (j3J)] can be separated from the spin operators. 

To compare with the measured line shifts, we calculate perturbatively the first moment, 
(H — H Ies ), where H Tes is the resonant field without the spin-spin interactions, and subtract 
from it the first moment at a reference temperature, 30 K. This reference temperature 
was chosen because at higher temperatures additional line shifts can be expected from the 
temperature dependence of D. The calculated line shifts to zero order in /k B T contain 
the exchange coupling constant J and the effective dipole field A = J2j^kAjk/N (N is the 
number of molecules in the sample and the summation runs over all molecules) as variable 
parameters. This zero-order result in Tv- 1 ' /k B T depends on J and A through the terms 
(S n i — £ n ) in Eq. (JJJ). Since A depends on the sample shape due to the field-induced net 
magnetization, Ewald's method 39 is not sufficient to estimate its value in our case. 

To compare with the measured linewidths, we need to calculate the second central mo- 
ment, ((H-(H)) 2 ), which is equivalent to ((H - H ves ) 2 ) - ((H - H ies )) 2 , where ((H-H ves ) 2 ) 
and ({H — H tes )) 2 are calculated perturbatively to zero order in /k B T. 37 The square root 
of the second central moment is proportional to the broadening due to the spin-spin inter- 
actions. The quantity ((H — H ves ) 2 ) includes the following six terms, Yfij Jfj> Yfij -hj\j-, 
Ylij A 2 j, J2'i jk JijJjk, J2i jk JijAjk, and J2' ijk AjA jk . The quantity ({H - H res )) 2 includes three 
terms, (£y Jij) 2 /N, ^OCEy Aj)/N, and {E'ij Aj) 2 / N . Here and Z' ijk run over 
all molecules with the constraint that no two indices in the summations may be the same, 
and Jij and Jj k are nonzero for nearest-neighbor molecules only. The coefficients of the last 




three terms (the summations over three indices i, j, and k) in ((H — H res ) 2 ) are the same 
as those of the three terms in ((H — H Tes )) 2 . To simplify the second central moment, we use 
the identities 

2 

= E'4+E%^, ( 9 ) 

ij ijk 

= E 'JijAij + 2J 'JijAjk , (10) 

ij ijk 

= E'4 + E'^A*- (ii) 

ij ijfc 

Summations over four different indices do not appear in Eq. (fTTj) because J2'ij Ay = 
N A\j by translational invariance. We thus have only three undetermined terms, 
H'ijJfji Y^'ijAfj, and J2 V - JijAij in the calculated second central moment. The exchange 
coupling constant J could be determined from the line-shift analysis. The remaining two 
terms, J2'ij A 2 j and JijAij, can, in principle, be calculated from the exact geometry of the 
system. However, in our study, we take the two terms as variable parameters, and compare 
their optimum values with the calculated values. Physical justification for this is provided 
in Sec. IV.A.2. 

We fix the EPR frequency at v — 116.9 GHz and vary the temperature from 4 K to 
35 K for the linewidth analysis (2 K to 30 K for the line-shift analysis). We do not analyze 
the experimental data above 35 K because at higher temperatures excited states (effective 
spin S < 10) might play a role. Only good-quality (high signal-to-noise ratio) experimental 
data were selected. For the line-shift analysis, the exchange constant J and effective dipole 
field A are varied, while for the linewidth analysis, T = J2'ijA 2 j/N, A = J2„- JijAij/ 'N, and 
the standard deviation of D are varied within experimentally acceptable ranges in order to 



fit the experimental data. Note that in Ref. 



20 the molecules of Fe« were assumed to be 



distributed on a simple cubic lattice in a spherical sample. For spherical samples the value 
of A is zero for all effective dipole distances. Thus only one parameter, either V or the 
effective dipole distance, sufficed in the linewidth analysis. The real samples, however, were 
not spherical. The Fes sample examined^ was a thin rhombic platelet with acute angle of 
about 60 degrees, edges of length 0.7 mm, and thickness 0.17 mm, and the molecules in 
Fe§ are distributed on a triclinic lattice. Therefore, we here use two fitting parameters (A 
and T) to consider the dipole-dipole distributions in the real experimental samples described 
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above. 



B. Mni2 

For M1112 we consider an effective single-spin Hamiltonian which satisfies tetragonal sym- 
metry, 

H = -DS 2 Z - CS A Z - 9HbH z S z (12) 

with D=0.55kB, C=1.17 x 10~ 3 A;b, and g=1.94.— Here we consider the case in which the 
applied field is along the easy axis (the z axis), and we neglect the small transverse fourth- 
order anisotropy term, S" 4 + S*. 

The technique is the same as for Feg, except for the following: (i) The ^-strain effect 
provides a weak temperature dependence to the linewidths, caused by the temperature de- 
pendence of the natural linewidths. (ii) The resonant field without the spin-spin interactions 
is modified to 

_ hv-D(2M s -l) + C(4M*-6M 2 +AM s -l) 

tl xes = • (i-O) 

(iii) To calculate the natural linewidths, we use the stren gth of the coupling between the 
spin system and a surrounding phonon heat bath from Refj4fJ. (iv) Dipoles are distributed 
on a centered tetragonal lattice with sample dimensions 1 x 0.15 x 0.05 mm 3 . (In Ref. [^j, 
the dipoles were assumed to be distributed on a body-centered cubic lattice, and the sample 
was assumed to be spherical.) (v) The easy anisotropy axis is along the long side of the 
needle-shaped sample, (vi) For Mni2 the measured line shifts are negligible compared to 
the measured linewidths, so that we do not have to consider the exchange interaction and 
the effective dipole field (J = and A = 0, so A = 0). Thus, the second central moment 
((H — (H)) 2 ), which is proportional to the measured linewidths, is identical to the second 
moment {{H — H res ) 2 ). 

To compare with the experimental data, the frequency is fixed at v — 189.123 GHz, 
and the temperature is varied from 10 K to 40 K. Our analysis ends at 40 K because at 
higher temperatures excited states (effective spin S < 10) might play a role . 41 ' 42 ' 43 Only 
EPR spectra of good quality were selected, and T and the standard deviations of D and g 
were varied within acceptable ranges in order to fit the experimental data. 
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IV. RESULTS AND DISCUSSION 



We show that the spin-spin interactions (dipolar and/or exchange interactions) alone de- 
termine the trend of the temperature dependence of the line shifts and the linewidths. From 
the line-shift analysis, we can estimate the orders of magnitude of the exchange interaction 
and the effective dipole field and obtain their signs. Using this information, we can also ex- 
plain quantitatively the linewidths, including the D-strain and/or the g-strain effects which 
give rise to a strong M s dependence but weak temperature dependence of the linewidths. 
The spin-spin interactions contribute more to the linewidths for Fe§ than for Mn 12 , mainly 
because the D-strain effect is dominant over the spin-spin interactions for Mn 12 , while it is 
comparable with the spin-spin interactions for Fe§. This explains the different temperature 
behavior of the linewidths for Fes and Mni 2 . The set of parameter values which best explains 
the experimental data has some systematic theoretical uncertainties, which are difficult to 
calculate exactly. 



A. Fe 8 

1. line shifts 

It is known that the value of the uniaxial anisotropy parameter D may vary smoothly 
with temperature.— To gauge the importance of this effect, we first assume that D has 
a temperature dependence such as Fig. [Ha). Then the line shift, (H)(T) — {H)(T =30 K), 
becomes monotonically temperature dependent, as shown in Fig.^b). Comparing Fig.^b) 
with the experimental data in Fig.|2fa), we see that the monotonic temperature dependence 
of D cannot by itself explain the complicated temperature dependence of the measured line 
shifts. It could presumably contribute together with other factors described below. However, 
we hereafter take Dasa temperature independent parameter, for the sake of simplicity and 
because the exact temperature behavior is not yet known. We also note that the distribution 
in D does not change the line positions. 

Next, we consider the effect of the spin-spin interactions between molecules on the mea- 
sured line shifts. If we first ignore the exchange interaction and consider the dipolar interac- 
tion only for a spherical sample with dipoles distributed on a simple cubic lattice (which is 
the assumption made in Ref 2(1) . then the effective dipole field A vanishes, so that there is 
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no line shift to zero order in TC^'/ksT. Higher-order corrections [the second term in Eq. (jSJ)] 
provide a much smaller and qualitatively different temperature dependence from that seen 
in the measurements [compare Fig. 0(a) with Fig. Ufa)]. If we include a non-zero effective 
dipole field A only in the zero-order calculation, then a negative effective dipole field moves 
the line shifts for all the transitions down below zero [Fig. E^b)]. On the other hand, if we 
include an exchange interaction only in the zero-order calculation, then a ferromagnetic ex- 
change interaction (negative J) moves the line shifts up above zero for all transitions except 
M s = +10 [Fig. 01(c)] . In both cases, the calculated line shifts behave very differently from 
the measured shifts [Fig. |2fa)]. Therefore, we need to include both the effective dipole field 
and the exchange interaction in order to explain the measured line shifts. An inter-molecular 
exchange interaction was recently observed for the different types of single-molecule mag- 
nets, Mn4 and Mn4 dimerj 4 ^ Since the effective dipole field A depends on the sample shape 
due to the field-induced net magnetization, we do not use Ewald's method^ to estimate A 
and thus leave it as a fitting parameter. For the exchange interaction, we assume that the 
coupling constant Jij is isotropic along the a, b, and c directions of the triclinic unit cell 
(although in experimental samples the exchange interactions are highly anisotropic), so that 
the coordination number is 6. The optimum values of A and J are A opt ~ —20 gauss and 
J op t ~ —7 gauss (~ 1 niK) respectively. With these optimum values, the calculated line 
shifts [Fig.Efb)] reproduce well the trends of the temperature dependence of the experimen- 
tal data [Fig. 121(a)] . Figure |2Kc) shows a direct comparison between theory and experiment 
for a few transitions. 

The negative sign of the effective dipole field (A < 0) indicates that dipoles are anti- 
ferromagnetically coupled. This result seems to be in conflict with the prediction that the 
dipolar Ising spin system with the same structure as Fes is ferromagnetically ordered . 4 ^ 4 ? 
However, as pointed out in Ref. |42l, the energy difference between the ferromagnetic and 
antiferromagnetic states is so small that any neglected effects may shift the ground state to 
an antiferromagnetic state. The negative sign of the exchange coupling constant (J < 0) 
corresponds to ferromagnetic coupling between the effective spins of the molecules. Thus, 
if there exists any ordering for Fes, then the ordering temperature should be estimated by 
considering both the exchange and dipolar interactions. The two interactions compete with 
each other, thereby reducing the possible ordering temperature to a lower value than the 
ordering temperature with only one of the two interactions considered. 
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Finally, we show the calculated line shifts with several other parameter values that are 
different from the optimum ones. If the effective dipole field and the exchange interaction 
both change signs, then the calculated line shift also changes its sign. If the sign of A is 
opposite to the sign of J, then the magnitude of the calculated line shift for M s = +10 is 
much smaller than those for the other transitions, which does not agree with the experimental 
data [compare Fig. OJd) with Fig. Ufa)]. Figure Efe) shows the calculated line shifts with 
J < J opt and A = A opt . Figure Off) shows the calculated line shifts with J = J opt and 
A < A opt . All three figures [Figs. E(d) -(f)] are significantly different from Fig. Efb) with the 
optimized values. 

2. linewidths 

For Fes the distribution in D and the spin-spin interactions contribute approximately 
equally to the inhomogeneous line broadening. Figure |U shows the calculated line broaden- 
ing due to the D-strain effect only as a function of temperature at v = 116.9 GHz. Here the 
standard deviation of the Gaussian distribution in D, <td, is approximately 0.0064D. The 
line broadening caused by the D-strain only becomes temperature dependent, because the 
natural linewidths depend on temperature. The distribution in D makes each molecule sub- 
ject to a slightly different resonant field. A measured line shape is a sum of many Lorentzian 
line shapes with a natural linewidth and different resonant fields. We can calculate the vari- 
ance of the resonant field, <t^(2M s — l)/gfi B , due to the distribution in D from the expression 
for H Tes . If the natural linewidths are comparable with the variance of the resonant field, 
then the effect of temperature is significant. If the natural linewidths are much smaller than 
the variance, then the effect of temperature is negligible. The natural linewidth at 10 K 
(35 K) varies from 7 G (29 G) to 79 G (235 G) as M s changes from +10 to +3, using the 
parameter values in Refl36l The variance of the resonant field, <7d(2M s — l)/gfiB, varies 
from 260 G to 70 G as M s is varied from +10 to +3 with a D = 0.0064D. Thus, we find 
that for small M s the natural linewidths are comparable with the variance, while for large 
M s the natural linewidths are much smaller than the variance. Therefore, for small M s the 
calculated linewidths show a substantial temperature dependence, while for large M s there 
is only a very weak temperature dependence (see Fig. HJ). 

In Fig. El the calculated line broadening caused solely by the spin-spin interactions at 
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fixed D is shown vs temperature at v = 116.9 GHz. Here we use the exchange con- 
stant, J = —7 gauss, which was estimated from the measured line shifts (Sec. IV A.l), 
r = J2'ijA^j/N = 86 gauss 2 , and A = Yfij JijAij/N = —156 gauss 2 . For the ground state 
M s = +10, the linewidths decrease with increasing temperature in the whole examined 
temperature range. For M s = +9, +8, and +7, the widths first increase sharply with tem- 
perature at low temperatures, and then decrease slowly with temperature at high temper- 
atures. For M s = +6, +5, +4, and +3, the widths increase with increasing temperature in 
the whole examined temperature range. As the temperature increases, the M s dependence 
of the line broadening due to the spin-spin interactions decreases. This trend was also seen 
in the experimental linewidths (shown as symbols in Fig. |UJ), confirming that the spin-spin 
interactions are essential to understanding the temperature dependence of the linewidths. 

The trends of the temperature dependence can be qualitatively understood through the 
relative magnitude difference of the thermal energy and the Zeeman energy splitting between 
the states M s = +10 and M s = —10. If the Zeeman energy splitting is much larger than 
the thermal energy (this occurs at low temperatures), then the system is polarized. Thus, 
higher temperature provides larger populations in higher energy levels within the same 
potential well where the ground state is located. This leads to an increase in the randomness 
of the spin orientation so that linewidths become larger with increasing temperature. If 
the Zeeman energy splitting is much smaller than the thermal energy (this occurs at high 
temperatures), then some energy levels in both potential wells are populated. In this case, 
thermal fluctuations increase rapidly with increasing temperature, so that the duration 
time of the local magnetic field due to neighboring molecules becomes shorter than the 
spin-spin relaxation time T 2 (the inverse of the natural linewidths). Eventually, at very 
high temperatures the local field is averaged out. Therefore, the linewidths decrease with 
increasing temperature, which usually occurs in paramagnetic materials with very small 
or zero single-ion anisotropy. 48 (This effect is called motional narrowing.^) Therefore, the 
"crossover" temperature where the maximum of the linewidth occurs must be proportional 
to the Zeeman energy splitting between M s = +10 and M s = —10, which is 2g/i£if res • 10. 
For example, for the transition M s = +10 — * +9 at v = 116.9 GHz, the resonant field is less 
than 0.1 T, so the Zeeman splitting is about Iks- Consequently, its crossover temperature is 
below the examined temperature range. The crossover temperature increases with decreasing 
M s because the resonances are observed at increasing fields for decreasing M s . For the 
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transitions M s = +9 — > +8, +8 — > +7, and +7 — > +6, the crossover temperatures are within 
the examined temperature range (see the inset in Fig. 6 of Ref J2 ll) . For M s = +6, +5, +4, 
and +3, the crossover temperatures are above the studied temperature range. 

Figure El shows the experimental data (symbols), and our calculated linewidths (curves), 
including both the D-strain effect and the spin-spin interactions with gd pa 0.0064Z), J = —7 
gauss, T = 86 gauss 2 , and A = —156 gauss 2 . (The spread in D here is different from 



that reported in Ref 



20l because gd is sample dependent and the samples examined were 



different. The value of V corresponding to an effective dipole distance of 12 A (Fe§) in 
Ref was about 203 gauss 2 .) As shown in Fig. our calculated linewidths agree well with 
the experimental data, except in the low-temperature range for large M s transitions (M s = 
+ 10, +9, +8). The experimental linewidths for M s = +10 show 1/T dependence in the whole 
examined temperature range.— For M s = +10, +9, and +8, the calculated linewidths are 
appreciably smaller than the experimental linewidths below 10 K. As a possible explanation 
for this discrepancy, we speculate that at low temperatures and large M s , (i) our assumption, 
TiS 1 ' /ksT "C 1, may break down, and/or (ii) there might be other linewidth contributions 
that we have neglected, which should be included along with the dipolar and exchange 
interactions. In principle, when Ti^/ksT is not much smaller than unity, a first-order 
calculation in TiS 1 ' /ksT produces corrections of 0(1/ ksT). But its implementation is quite 
complicated, and a first-order calculation may anyway not be sufficient to explain fully the 
measured linewidths. 

Introducing the concept of the crossover temperature to explain the temperature de- 
pendence of the widths seems to be successful for v = 116.9 GHz. However, recent EPR 
experiments (Fig. 8 in Ref 21) showed that even when the frequency increased to v — 145.9 
GHz such that the resonant field for the ground state transition (M s = 10) is approximately 
1 tesla, the linewidths for this transition still increased with decreasing temperatures down 
to 2 K. This cannot be explained using the reasoning given above, because the crossover 
temperature for the transition is approximately 15 K, so that the linewidths should decrease 
with decreasing temperature below about 15 K. At present, we do not fully understand the 
broadening of this ground-state transition. (At v = 145.9 GHz, for other transitions than the 
ground-state transition, the temperature dependence of the linewidths can be understood 
using the concept of the crossover temperature.) 

Consideration of exchange interaction in the linewidths slightly reduces gd (from 0.0076D 
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to 0.0064D) and the dipolar interaction (from T = 103 gauss 2 to T = 86 gauss 2 ). However, 
the quality of the linewidth fit including exchange interaction is comparable to that without 
exchange interaction since the exchange coupling constant is very small compared with the 
linewidths.— The two fitting parameters, T = 86 gauss 2 , and A = —156 gauss 2 , can be 
calculated using the exact geometry of the system. The calculated values are T ca ; = 500 
gauss 2 and A ca ; = —137 gauss 2 , when the easy axis is 9° off from the a axis toward the 
positive b axis, and 7° off from the ab plane.— The optimum value of A is quite close to 
A ca ;, in contrast with T. Possible reasons that the optimum value of T is much smaller 
than r ca ; are as follows: (i) In our calculation, we considered each molecule to be a point 
dipole. If we consider the atomic positions of the eight Fe ions in each molecule and calculate 
the dipolar interaction between Fe ions in different molecules, then the sum of the squared 
dipolar interaction, T, can be significantly reduced, (ii) Recent NMR experiments for the 
single-molecule magnet Mn 12 showed some spin-density leakage onto the ligands . 27 i 28 i 29 This 
indicates indirectly that for Fes the spin density in a single molecule may not be confined 
only on the core, which would thus reduce the magnetic moments of the eight Fe ions. The 
above two reasons are our speculations to explain the discrepancy, but it is still unclear 
why considering the atomic structure within each molecule does not substantially change 
the value of A. 

B. Mni2 

The experimental data for Mn 12 are limited to resonance frequencies below 190 GHz, so 
that the large M s transitions where the line shifts are significant cannot be observed for 
these low-frequency measurements. Additionally, the linewidths for Mn 12 are an order of 
magnitude larger than those for Feg. Therefore, relatively small line shifts are probably 
masked for Mni2. Thus, hereafter, the small line shifts are ignored in our analysis, so the 
exchange interaction and the effective dipole field A need not be considered in the linewidth 
analysis for Mn 12 . The sources of the line broadening are then D-strain, g-strain, and the 
dipolar interaction. 

Line broadening caused by the D-strain and ^-strain effects for Mni 2 is found to have a 
weak temperature dependence (not shown), which is similar to the line broadening due to 
the D-strain effect for Fes- The contribution of the dipolar interaction to the linewidths is 
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shown vs temperature in Fig. [7| Here T = J2'ij Af^/N = 203 gauss 2 . The dipolar broadening 
increases with increasing temperature for M s = +6, +5, +4, +3, and +2. We do not see the 
regime where the dipolar broadening decreases with increasing temperature, because the 
crossover temperature for M s = +6, about 32 K (the resonant field is about 1.6 T), is close 
to the highest temperature analyzed (40 K). Unlike Feg, the M s dependence of the dipolar 
broadening does not decrease with increasing temperature (the curves are almost parallel). 
This is also observed in the experimental data (shown as symbols in Fig. ISJ). 

We combine the three effects (.D-strain, g-strain, and dipolar interactions) to find that 
the calculated linewidths agree well with the measured linewidths with <jd ph 0.018D, o g pa 
0.002#, and T = 203 shown in Fig. |HJ (The value of T corresponding to an 

effective dipole distance of 14 A for Mn 12 in Ref|20j turned out to be the same as that 
obtained for Fes i n Ref. |20j) Here the standard deviation of g is quite small, so that we 



cannot rule out the possibility of a g = 0. Note that op and a g vary from sample to sample. 
The optimum parameter values found here are different from those estimated in Reffl 
because the examined Mni2 sample was different. The calculated value for T, with the 
exact geometry of Mni2 from RefQj (with each molecule considered as a point dipole), is 
F C ai = 397 gauss 2 ; this is, again, quite a bit higher than the optimum value for T, probably 
for the same reasons as for Fes although the origin of this discrepancy remains unclear. 
Overall, the temperature dependence of the linewidths for Mn 12 is weaker than for Fe§, 
because the distribution in D for Mni2 is roughly three times as wide as for Fes, and the 
dipolar broadening for Fes is comparable to that for Mni2. Thus, the distribution in D 
conceals the significant temperature dependence of the dipolar broadening for Mni2. 

As a consistency check, we also used the same values of the three parameters (op ~ 
0.018.D, o g pa 0.002(7, and T = 203 gauss 2 ) to analyze the measured linewidths2i as functions 
of the energy level M s for several frequencies [y = 127.8, 148.5, 169, 181.8, and 189.1 GHz) 
at a fixed temperature (T=20 K). Our calculated linewidths are in good agreement with the 
experimental data as shown in Fig. El Due to a dominant contribution of the distribution 
in D to the linewidths, the linewidths do not depend much on the resonance frequency. 
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V. CONCLUSION 



We have investigated how the EPR line shifts and linewidths vary with temperature 
for different energy levels M s with the applied field along the easy axis for the single- 
molecule magnets Fe§ and Mn 12 . Our calculations consider the spin-spin interactions be- 
tween molecules, as well as distributions in D and g. We have found that the distributions 
in D and g provide a weak temperature dependence to the linewidths, and that the spin- 
spin interactions (exchange and dipolar interactions) dominate the temperature dependence 
of the line shifts and the linewidths. For Feg, the line-shift analysis (Figs. El and EJ) pro- 
vides possible evidence of an exchange interaction between molecules, and it determines the 
sign and order of magnitude of the exchange interaction. The competition of the suggested 
exchange interaction with the dipolar interaction would tend to lower a possible magnetic 
ordering temperature. A small exchange interaction does not affect the linewidth analy- 
sis significantly because the exchange coupling constant is much smaller than the typical 
linewidths. Table H] summarizes the optimized values of the parameters used in our analysis. 
Those parameters are, in principle, independent of the resonance frequencies, but some of 
them are expected to be somewhat batch and shape dependent, as indicated in Table |TJ 
Because of the much broader distribution in D for M1112, the linewidths for Fes show a 
stronger temperature dependence than those for Mn 12 . This conclusion also corroborates 
our assumption that D is distributed for both materials,— although the microscopic origin 
of this spread is not yet well understood . 23124 
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TABLE I: Optimum values of the parameters used in the line-shift and linewidth analysis. Here 
od is the standard deviation of D, a g is the standard deviation of g, J is the exchange coupling 
constant between nearest-neighbor molecules (negative sign means ferromagnetic interaction), A = 
J2ijAij/N which vanishes for spherical samples, T = J2ijAjj/N, and A = J2ijJijAij/N. With 
experimentally determined values of D and gr^ f° r F e 8 we optimize J and A for the line shifts 
and ctd, T, and A for the linewidths, while for Mni2 we optimize ao, &g, and T for the linewidths. 
The parameters are essentially independent of the measurement frequencies, but some of them are 
expected to be somewhat batch and shape dependent, as marked by X in the table below. 



Fe s 



Mn 



12 



Batch 



Shape Size Crystal structure 



D 

g 

J 
A 

r 

A 



0.288 k B 
0.0064D 
2.00 

—7 gauss 
—20 gauss 
86 gauss 2 
— 156 gauss 2 



0.55 k B 
0.018D 
1.94 
0.002 5 



X 

X 
X 



203 gauss 2 



X X 
X X 
X X 



X 
X 
X 



"Batch, shape, and crystal structure independent 
h Batch dependent 

c Shape, size, and crystal structure dependent 
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FIG. 1: (a) Hypothetical smooth temperature dependence of the uniaxial anisotropy parameter D 
for Fe 8 . The functional form used here is D(T) = -0.710665 + exp{-l/[10(50 - T) + 250]}. (b) 
The resulting calculated line shifts (the peak position at a given temperature minus the position 
at T = 30 K) due to this temperature dependence of D, shown vs temperature at v = 116.9 GHz. 
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FIG. 2: (a) Measured line shifts vs temperature at v = 116.9 GHz for Fes- (b) Calculated line 
shifts vs temperature at v = 116.9 GHz for Fes- Here the calculation is performed to zero order in 
H<ti/k B T. We use the effective dipole field A = ^ Aj/N = -20 gauss and the exchange coupling 
constant J = —7 gauss, (c) Measured line shifts (symbols in (a)), superimposed on calculated line 
shifts (curves in (b)), for several transitions for comparison. 
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FIG. 3: Calculated line shifts vs temperature at v = 116.9 GHz for Fes, with (a) only the dipolar 
interaction for a spherical sample in higher-order calculations ( J = 0, A = Yl'ij Aij/N = 0, and 
r = T,ij A ij/ N + 0) (b) only the effective dipole field, A = -20 gauss, (c) only the exchange 
coupling constant, J = —7 gauss, (d) A = —20 gauss < and J = +7 gauss > 0, (e) A = —20 
gauss and J = — 12 gauss, and (f) A = —30 gaMs and J = —7 gauss. 
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FIG. 4: Calculated full width at half maximum (FWHM) caused by the Gaussian distribution in 
the uniaxial anisotropy parameter D only, shown vs temperature from 4 K to 35 K at v = 116.9 
GHz for Fe§. M s = 10 indicates the transition from the energy level M s = 10 to M s = 9, etc. The 
standard deviation of D is approximately 0.0064D. 
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FIG. 5: Calculated FWHM caused by the spin-spin interactions only, shown vs temperature at 
v = 116.9 GHz for Fes- Here the exchange constant J is —7 gauss, V = J2ijAfj/N = 86 gauss 2 
and A = JijAij/N = -156 gauss 2 . 
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FIG. 6: Calculated (curves) and measured FWHM (symbols) vs temperature at v = 116.9 GHz 
for Fes- Here we use the standard deviation of D, od ps 0.0064D, the exchange coupling constant, 
J = — 7 gauss, r = 86 gauss 2 , and A = —156 gauss 2 . The solid curves, from bottom to top, 
correspond to M s = 3, 4, 9, 10. See the text for possible sources of the discrepancy between 
theory and experiment for low temperatures and large M s . 
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FIG. 7: Calculated FWHM for Mni2 caused by the dipolar interactions only, shown vs temperature, 
at v = 189.123 GHz, with the sum of the squared dipolar interactions T = 203 gauss 2 . The 
examined temperature range for Mni2 is from 10 K to 40 K. 



26 



2000 



1000 



Mn 12 










Ms=6 


[] 








n 


□ 








u 








ii 


[] 

■ ■ 


□ 
■ 


■ 


■ ■ 

■ 


■ 


■ 




ii — 










Ms=5 

— — 1 ~ 




o 



• 


•e- 
• 


• 


• • * 


O 
• 


Ms=4 
• 

Ms=3 




o 


o 












Ms=2 





10 20 30 40 

Temperature (kelvin) 



FIG. 8: Calculated (curves) and measured (symbols) FWHM vs temperature at v = 189.123 GHz 
for Mni2- Here the D-strain (ctd ~ 0.018-D), ^-strain {a g ~ 0.002^), and the dipolar interactions 
(r ss 203 gauss 2 ) are included in the calculated linewidths. 




FIG. 9: Calculated (curves) and measured (symbols) FWHM vs energy level M s at T = 20 K for 
v = 127.8, 169, and 189.1 GHz for Mni2- Here the values of <td, a g , and T are the same as those 
in Fig. |HJ Because of a relatively small contribution of the dipolar interaction to the linewidths, 
the linewidths do not change much with the resonance frequency. 
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